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Abstract 

Wave scattering from rough surfaces in addition the inverse scattering is an interesting approach 
to obtain the surface topography properties in various fields. Analytical expression in wave scattering 
from some known rough surfaces, not only help us to understand the scattering phenomena, but also 
would prove adequate to be a criterion to measure the information for empirical rough surfaces. For 
a rough surface with an exponential height correlation function, we derive an analytical expression 
for the diffused part and expanded it in two asymptotic regimes. We consider one surface as slightly 
rough and the other as very rough based on the framework of the Kirchhoff theory. In the end, we 
have measured the role of various Hurst exponents and correlation lengths on scattering intensity in 
self-affine surfaces. We have shown that by increasing the Hurst exponent from H — to H = 1, 
the diffuse scattering decreases with the scattering angle. 

PACS: 42.25.Fx, 68.35.Ct 

1 Introduction 



The study of rough surfaces is needed in many scientific fields. This is due to the fact that roughness is 
a scaling parameter and depending on the scale of observation which plays an important role in natural 
phenomena. The behavior of scattered fields from rough surfaces is of interest to may research topics, 
ranging from X-ray to radar wave in a long time in various fields. One of the most general approaches 
to study wave scattering is the Kirchhoff theory [I] [21 (3J [H [3] . This theory is an electromagnetic theory 
and is known as a tangent plane theory which is most widely used to calculate the distribution of the 
specular and diffuse parts of the reflected light. The Kirchhoff theory treats any point on a scattering 
surface as a part of an infinite plane, parallel to the local surface tangent [3] . 
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There are two types of problems in this area: (a) direct problem, (b) inverse problem. The direct 
problem is concerned with determining the scattered field from the knowledge of the incident field and 
the scattering obstacle [6j [7j |8j [9j [10l [IT]. Interesting works in the context of wave scattering has 
been carried out Brewster's scattering angle [12] and the Rayleigh hypothesis [13] have been studied. In 
addition, Ingve et al P3] studied wave scattering from self-affine surfaces where Leskova et al [JJ>] studied 
the coherence of p-polarized light scattered from rough surfaces. Since the wave-length could be either 
smaller or greater than the height fluctuations, the two scale theory was proposed [16 . In a further 
study, effects of interference of two beam scattering was considered [17] . Scattering from multilayer is 
studied by Carniglia [18]. Elson et al [19] investigated vector wave scattering theory to estimate the 
angular distribution of scattered light from optical surfaces. Some of the researcher studied shadowing 
effect on rough surfaces [201 121] . The inverse problem is concerned with inverse scattering techniques 
to [551 HS1 [23 US] measure the statistical properties of rough surfaces. Sinha et al showed how scaling 
properties exists in X-rays and neutrons from rough surfaces [35], Jafari et al [37] obtained the surface 
roughness, Dashtar et al used inverse wave scattering to determine the height distribution on a rough 
surface [28] . 

The results indicate that when the wavelength is larger than the correlation length, the scattered 
wave differs with the case where the wavelength is smaller than the correlation length. In this work an 
analytical expression for the diffused intensity is obtained. For the case where the height correlation 
function has an exponential behavior. In addition scattering from self-affine surfaces is studied for various 
roughness exponents and correlation lengths. 

The paper is organized as follows. In Sec. II we present a theoretical description of Kirchhoff theory 
and surface roughness is the source of scattered field. The self-affine fractal model of the surfaces in Sec. 
III. The comparison of coherent and diffuse parts of scattered intensity and scattering from self-affine 
surfaces with different parameters are described respectively in Sec. IV. Finally, some general conclusions 
are presented in Sec. V. 
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Figure 1: Schematic figure to show the geometry used for wave scattering from a rough surface. 

2 Kirchhoff theory 



In Kirchhoff theory, the incident held is written as ip mc (r) = exp(— iki nc .r), where k and r are the wave 
number and position respectively. Fig. 1 shows schematically the scattering phenomena from a rough 
surface with dimensions —X < xq < X, —Y < ya < Y. The scattered field in r is denoted by ip sc 
and relies on three important assumptions: a) The surface is observed in the far field, b) No point on 
the surface has infinite gradient. Therefore, the Fresnel laws can be locally applied, c) The reflection 
coefficient, Ro, is independent of the position on the rough surface. The total scattered field over the 
mean reference plane Am, is given by [3], 

ike ikr f dh dh 

V c {r) = —. / [a- lb- c) x exp(ik[Ax + By Q + Ch(x ,y )})dx dyo, (1) 
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c = cos6> 2 (l + Rq) - cos6»i(l - R Q ). 
The coherent and the diffuse intensities are given by [3] , 

Icoh = he~ a , 

L2 p2 poo 

<I d> = 7T^ AMe ~ a / Jo{kR\/ A 2 + B 2 ) x [ e s Cor W - l]RdR, (2) 
^ r Jo 

where Iq is the scattered intensity from a fiat interface. The parameter g is equal to = k 2 u 2 C 2 , where 
C = cos#! + cos #2, and the height- height correlation function Cor(R) which is definite for an isotropic 
surface ((h(x)) = 0) is equal to W*+W)> 

3 Scattering from self-affine surfaces 

3.1 Self-affine surface 

The correlation function of self-affine surfaces is known by correlation length, £, which is the mean 
lateral length of surface features, e.g. grain size or other. Following Sinha et al. ES GU [30], can be 
approximated by an analytic function, Cor(R) w e v £ ' , which contains the following limits. For 
R ^> £, the correlation vanishes, (Cor(R) = 0) and for R <C £, the correlation function would look like 
Cor(R) ~ 1 — (■^) 2 ' ff , and roughness exponent H would be between zero and one. The large values of 
H correspond to smoother height-height fluctuations, while small values of H characterize irregularity 
in height for rough surfaces at the length scales (R <C £) [2§1 [3D] ■ 

For a Gaussian correlation function (H = 1), Eq. (2) yielded an expression for the diffuse scattered 
intensity [3]: 

k 2 F 2 ee-v k 2 (A 2 + B 2 )e , ... 

n— 1 

3.2 Analytical expression for a rough surface with exponential height corre- 
lated 

Knowing any exact results in wave scattering from rough surfaces is helpful to better understand the 
scattering phenomena. One of these surfaces is exponential height correlated surface which could be an 
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Figure 2: (Color online) Dependence of diffuse scattered intensity on ka for different correlation length 
(£ = 500, 1000, 1500nm) and for angles 9i = 9 3 = 0°, 2 = 20°, wavelength A = 500nm and H = 0.5. 

important and suitable criteria to measure information content in other rough surfaces. We obtain an 
analytical expression for the correlation function Cor(R) — e ( i ' with H = 0.5, which appears in an 
exponential form Cor(R) = e~ T. With regarding e 9Cor ^ — J2^=o 9 c ° r n \ ^ , the diffuse part of the 



scattered intensity will be: 



< Id > = 



A M e- g / J a {kR^/A 2 + B 2 ) x ^ 



n=l 
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(4) 



Eq. (4) enables the study of rough surfaces in two asymptotic limits: 

A. Slightly rough surfaces, For slightly rough surfaces (g <C 1), the series in Eq.(4) converges 
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quickly and only the first term needs to be considered. So, the diffuse field intercity becomes: 

k 2 F 2 £ 2 

< Id >= ^-^A M e 9 g- — — 3 . (5) 

2irr 2 {l + k 2 (A 2 + B 2 )i 2 )2 

B. Very rough surfaces, The total intensity is (I) = (I co h) + (Id) = g 9 Iq + (Id), when the surface 
is very rough (g 3> 1), the coherent field will be negligible and the total intensity is equivalent to the 
diffuse ones [3] (Chap. 4, Page 92). In addition, we know that I oc X. , when sample size is 

larger than wavelength, the coherent field appears just in specular angle. For the non-specular angle we 
have only diffuse field (A < X,Y or kAX > 1 X. S "^ X) = n.S(kA)) [31]. So, according Eq.(2) the 
total scattered intensity which is approximately equal to one could be written as: 



2 tp2 



< I >= —^A M / J (kRVA 2 + B 2 )e-3(e 9CoriH > - l)RdR. (6) 
27rr^ J 



When the length scale is larger than the correlation length (R 3> £,), the correlation function vanishes 
(Cor(R 3> £) = 0). In this case, according to Eq. (6), there is no scattered intensity in the large length 
scale (R 3> £). For convergence, it is sufficient to extend the numerical integration from zero to 5£. In 
the category which g ^> 1 and for small range of R, Eq. (6) is equal to: 

k 2 F 2 



<I> = 



2nr 2 
k 2 F 2 



A AI / JoikR^A 2 + B 2 )e~ g[1 - Cor{R)] RdR 
Jo 

/>oo 

A M / J (kRy/A 2 + B 2 )e~ 9 f RdR 
Jo 



2nr 2 

k 2p2 /v 



2nr 2 AM (i + k'(A'+B')(' ^ g2 - ( ? ) 

4 Effects of the correlation length and Hurst exponent on scat- 
tering intensity 

4.1 Variation of Correlation Length 

According to Eq. (2) the correlation function only effects the diffuse part leaving the coherent part 
uncffccted. Fig. 2 compares the diffused scattering vs ka, with a fixed A = 500 nm, for three values of 
the correlation length f = 500, 1000, 1500 nm and the angles (0 X =63 = 0°, 6 2 = 20°). We have chosen 
the value of parameters in the way which are compatible with experimental setup. The numerical results 
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which were obtained from Eq. (2) are shown for ka = to 4. A characteristic feature seen in Fig. 2 is 
the presence of a maximum in the diffused part as a function of the ka, which depends on correlation 
length. By increasing the correlation length the peak of the diffused part is shifted to large ka. Indeed, 
small fluctuation have low effects on larger wavelengths. Note that decreasing wavelength and increasing 
rms have the same effect on wave scattering. Figuratively speaking, a surface with shorter correlation 
length appears rough to small wavelengths and smooth to large wavelengths. Therefore, in the context 
of scattering, ka is more suitable than k or a alone. In other words the longer angles would play a more 
significant role on diffused intensity. This result in a decrease\increase of intensity for smaller\longer 
angles. In addition, the smaller correlation length changes its behavior as it shifts to larger ka. 

4.2 Variation of Hurst exponent 

The Hurst exponent appears in the correlation function and the correlation function enters into the 
diffuse part. So, the variations of Hurst exponent could effect the diffuse scattered intensity, too. To 
more sense the role of Hurst exponent in wave scattering, we study diffuse scattering in two regimes. Fig. 
3 depicts the variation of the diffuse scattering intensity with respect to scattering angle for different 
Hurst exponents for the two cases: A < £ & A > £. The first point in Fig. 3 is that by increasing 
H, the diffuse scattering intensity is reduced for angles close to the specular angle widening the curve. 
Thus it can be observed in larger angles. For non-specular scattering angle far from specular angle, this 
behavior is opposite, so by increasing H, the diffuse scattering intensity is increased. Another point is 
that for A < £ the fluctuations of the surface are observed better and the roughness is increased so the 
diffuse scattering intensity increases (Fig. 3a). For the case A > £, smaller fluctuations of the surface 
are not observed and surface seems smoother and the diffuse scattering intensity decreases sharply. This 
behavior can be seen by comparing Figs. 3a and 3b. Variation of the diffuse scattering intensity with 
respect to H has the same behavior for both wavelengths smaller and larger than the correlation length. 

There are two kinds of behaviors in scattering intensity based on wavelength smaller or larger than 
the correlation length. When the wavelength is less than correlation length, the small height fluctuation 
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is being observed. So the wavelength acts as a scale for the observation of the surface. The more 
observed height fluctuation leads to the more diffuse scattered intensity. There is opposite behavior 
when the wavelength of incident wave is larger than the surface correlation length. In this case, the 
incident wavelength is not able to recognize small fluctuations on the surface and loses the scale of 
scattering angle. This means that the surface seams smooth for this scale of observation (wavelength). 
This results in a decrease in the diffuse intensity. 

Fig. 4 shows the variation of the diffuse scattering intensity with respect to ka for different H for 
the non-specular scattering angles (0\ = 9% = 0° and 8 2 = 20°). By increasing iJ, the diffuse scattering 
intensity decreases, which is in good agreement with Fig. 3. In Fig. 4, by increasing H, the curve shifts 
to larger ka. 

In order to investigate the role of the Hurst exponent, a surface with the observed values of for the 
correlation length (£ = lOOOnm), and wavelength (A = 500nm) is considered. When ka is small, either 
by means of small a or large A, the scattering is not sensitive to the corrugation, independently of the 
Hurst exponent. For the chosen correlation length of 1000 nm, the diffuse intensity in Fig. 4 is similar 
for ka < 0.4, the influence of the Hurst exponent becomes apparent only for large ka. 

5 Conclusion 

The roughness of a surface effects the scattered wave from rough surfaces. We have found an analytical 
expression for spectral scattered intensity for a rough surface with exponential height correlation function 
and investigated the results for asymptotic regimes of slightly rough and very rough surfaces. The exact 
solution helps us understand the scattering phenomena and could be a suitable criteria to measure the 
value of information in a unknown rough surface. 

The role of the Hurst exponent has been investigated for self-affine rough surfaces. Such an examina- 
tion was performed over a wide range of surface topographies, from logarithmic (H = 0) to a power-law 
self-affine rough surface, < H < 1. The roughness exponent H has a strong impact on the diffused 
part of wave scattering mainly for relatively large correlation lengths. Therefore, Hurst exponent must 
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be taken carefully into account before deducing the roughness correlation lengths from wave scattering 
measurements. 
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Figure 3: (Color online) Dependence of diffuse scattered intensity on O2 for angles 6\ = 0°,0s = 0° and 
different Hurst exponent, correlation length £ = lOOOnm, standard deviation a = 50nm and the incident 
wavelength (a) A = 500n?7j, (b) A = 5000nm. 



12 




Figure 4: (Color online) Dependence of diffuse scattered intensity on ka for angles 0\ = O3 = 0° and 
9-2 = 20° different Hurst exponent, correlation length £ = lOOOnm, wavelength A = 500nm. 
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